THE SELF-GRAVITATING GAS IN THE PRESENCE OF DARK ENERGY 
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The non-relativistic self-gravitating gas in thermal equilibrium in the presence of a positive cos- 
mological constant A (dark energy) is investigated. The dark energy introduces a force pushing 
outward all particles with strength proportional to their distance to the center of mass. We consider 
the statistical mechanics of the self-gravitating gas of N particles in a volume V at thermal equilib- 
rium in the presence of A. It is shown that the thermodynamic limit exists and is described by the 
mean field equations provided N, V — > oo with N/Vz fixed and AVs fixed. That is, A — > for 
N, V — > oo. The case of A fixed and N, V — > oo is solved too. We solve numerically the mean field 
equation for spherical symmetry obtaining an isothermal sphere for A > 0. The particle distribution 
turns to flatten compared with the A = case. Moreover, the particle density increases with the 
distance when the cosmological constant dominates. There is a bordering case with uniform density. 
The density contrast between the center and the boundary may be significatively reduced by the 
dark energy. In addition, the critical point associated to the collapse (Jeans') phase transition is 
pushed towards higher values of N/[T V"3] by the presence of A > 0. The nature and the behaviour 
near the critical points is not affected by the presence of A > 0. 
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I. INTRODUCTION 

The self-gravitating gas in thermal equilibrium has been thoroughly studied since many years 0, 0, H 0, As a 
consequence of the long range attractive Newton force, the selfgravitating gas admits a consistent thermodynamic 
limit N, V — > oo with -^j- fixed. In this limit, extensive thermodynamic quantities like energy, free energy, entropy 

are proportional to N\jify USQ. 

We investigate in this paper how the presence of a cosmological constant affects the properties of the nonrelativistic 
sclf-gravitating gas in thermal equilibrium. 

The cosmological constant modifies the Poisson equation adding a negative term in the r. h. s. This introduces 
an extra force pushing outwards all particles with a strength proportional to their distance to the center of mass. 
This outward force is the nonrelativistic version of the exponentially fast expansion in the de Sitter universe. It is an 
antigravitational effect. 

We consider the statistical mechanics of the self-gravitating gas in the presence of a cosmological constant and show 
that it is described by the mean field equation in the thermodynamic limit. We derive the mean field equation in the 
thermodynamic limit of large number N of particles and large volume V with N/V 1 / 3 and A7s fixed. Furthermore, 
we consider the case A fixed for N and V going to infinity with N/V 1 / 3 fixed. 

As known, the usual thermodynamic limit N, V — ► oo with N/V fixed leads to collapse in absence of the cosmological 
constant, while the dilute limit keeping N/V 1 / 3 fixed leads to a gaseous phase Q,y,|4|- Here, we consider a vanishing 
A for V — > oo keeping A V 3 fixed as well as N/V 1 / 3 . 

The meaning of the condition A V 3 for TV, V — > 00 goes as follows. The dark energy in a volumen V is of the order 
V A while the thermal (kinetic) energy (as well as the interaction energy between the particles) is of the order N T. 
Since N ~ V 1 / 3 , A fixed implies that the ratio of dark and thermal energies is fixed for N — > 00. 

When A and N/V 1 / 3 are kept fixed for N, V — > 00 the cosmological constant dominates over self-gravity. The 
result is a set of non-gravitating particles in a harmonic oscillator with negative squared frequencies. This system is 
exactly solvable as we show. 

It is convenient to define the dimcnsionless parameter: 

^2AGm|. (1.1) 

All physical magnitudes are expressed in terms of £ and 

r] = Gm 2 N/[TL] , 

the dimensionless parameter used in ref . [2j . 

We solve the mean field equations numerically for spherical symmetry. That is, we find the isothermal sphere in 
the presence of a non-zero cosmological constant. 

We find that this 'antigravitational effect' of A flattens the particle distributions in thermal equilibrium. The 
density contrast (ratio of the density at the center and at the surface of the sphere) decreases for increasing A. The 
particle density at the boundary is conformally translated as a function of r\ by the presence of A (see fig. 0). 

In absence of cosmological constant the particle density is a monotonically decreasing function of the radial distance 

P -0. We find in the presence of a nonzero A that the particle density may be either monotonically decreasing 
or monotonically increasing with R. The particle density grows with R when A dominates over the selfgravitational 
attraction. 

In the special case rj = £ the gravitational self- interaction exactly compensates the cosmological constant. As a 
result the gas becomes a perfect gas with uniform density. 

The position of the collapse (Jeans') phase transition is pushed by the cosmological constant to higher values of 
rj. Namely, lower temperatures or volumes or higher number of particles. The nature of the phase transition is not 
affected by the presence of A > and the behaviour of the physical magnitudes does not change far from the phase 
transition (see figs. 01 |S] and EJ) . Far from the collapse point, the energy and the free energy are basically shifted by 
the presence of A > 0. 

The thermodynamic quantities: free energy, energy and entropy turn to be proportional to the number of particles 
in the thermodynamic limit. We compute them as functions of rj and £. We also consider the limiting case £ 3> rj 
which is solved analytically in sec. V. This corresponds to keep A fixed for V — *■ 00. 

As for A = the gas collapses in the canonical ensemble when the compressibility diverges and becomes negative. 
This happens for a value of rj smaller than the one where the specific heat is singular. 

This article is organized as follows: in sec. 2 we recall the equations of motion for nonrelativistic particles in the 
presence of a cosmological constant. In section 3 we present the self-gravitating gas in the presence of A > 0; first in 
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a hydrostatic approach, then the statistical mechanics and finally the mean field approach including the spherically 
symmetric case with the calculation of the physical magnitudes (energy, free energy, entropy, surface density, pressure 
contrast, specific heat and compressibility). In sec. 4 we present the physical results for the the self-gravitating gas in 
the presence of A > 0: phase structure and particle density behaviour in the various cases. In sec. 5 we consider the 
exactly solvable case where the cosmological constant dominates overwhemly over the self-gravitation. We present 
our conclusions in sec. 6. 



II. NONRELATIVISTIC GRAVITY IN THE PRESENCE OF THE COSMOLOGICAL CONSTANT 

The energy-momentum tensor for massive particles plus dark energy (the cosmological constant) takes the form 

T| = A 81 + p S aQ Spa , (2.1) 

where A stands for the cosmological constant, p for the energy density of the massive particles, < a, [3 < 3 and we 
assume nonrelativistic matter so its pressure can be neglected compared with its rest mass. The Einstein equations 
take the form 

B% = ftoG(n- 5 -%T} , (2.2) 
where i?2 stands for the Ricci tensor and T = T£. We find from ea. l|2.1|l that 

T = 4A + p and T ° - ~ T = \p - A . 

We are interested in the non-relativistic limit where gravitational fields are weak and we have for the metric 

ffoo = 1 + 2 V , g ik = -&ik , 
where V stands for the gravitational potential. One has for the 00 component of the Ricci tensor^ 

R% = V 2 V . 

Therefore, the 00 component of the Einstein equations lj2.2|l becomes @ 

V 2 V = AnGp- 8vrGA . (2.3) 

For zero cosmological constant we recover the usual Poisson equation, as it must be. Eci. (|2.3|) can be written in an 
integral form, 

V® = V«M + VW0 - -G j 4?^f - ^ e , (2.4) 

J \q-<r\ 3 

where V se i f is the contribution of the massive particles and Vdark is the contribution of the dark energy. Here we 
choose the center of mass of the particles as origin of the coordinates. Linear terms in q of the type 

a-q (2.5) 

with a an arbitrary constant vector, are ignored in V(q) since they correspond to a constant external gravitational 
field unrelated to the cosmological constant. [Notice that such terms have zero laplacian and are not determined by 
ea. (|2.3() ]. Such external field will accelerate the particles in the direction a. Changing the origin of coordinates, terms 
linear in q can be added to ea. H2.4jl . Additional constant terms are clearly irrelevant. 

The mass density takes the following form for a distribution of point particles at rest, 

i 

Here, m, stands for the mass of the particle at the point ql. The gravitational potential thus becomes, 



irr-\ mi 4ttGA 
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and the gravitational field at the point q, 



J-VKfl-^^^l (2.6) 

The potential energy of such a system of non-relativistic self-gravitating particles in the presence of a cosmological 
constant A thus takes the form, 

i<j i 

Therefore, the Hamiltonian can be written as 



2 m 

where pi stands for the momentum of the i-th particle. The cosmological constant contribution to the potential 
energy decreases for increasing values of the particle distances ri to the center of mass. As stated above, in absence 
of external fields, the terms like ea. (|2.5I) are absent in the potential energy. Anyway, in the center of mass frame such 
linear terms identically vanish upon summing over the particles in ea. ()2.7JI . 

Therefore, the gravitational effect of the cosmological constant on particles amount to push them outwards. Equiv- 
alently, the last term of the gravitational field eq. (|2.6[) points outward. It can be then said that the cosmological 
constant has an antigravitational effect. 

For simplicity we shall consider all particles with the same mass m. The case of unequal masses may be worked 
out generalizing the work in ref.jj] to nonzero cosmological constant. 

Dark energy introduces a kind of centrifugal force always pushing out. Notice that the last term in ea- H2.7|) has the 
structure of a rotational energy but with negative sign. Also, it contains the distance squared of the particles to the 
origin and not the distance squared to an axis as it is the case for rotational energy. 

It is useful to consider the case of just one particle in the presence of the cosmological constant in order to obtain 
physical insights. We obtain from eas. H2.7|) and (|2.8|) . 

8vrGA _ 
1=^1- 

This is an oscillator equation with imaginary frequency and solution, 

q(t) = <f(0) cosh Ht + — <f(0) sinhfft , (2.9) 
H 



where H = y 87r ^ A ■ The particle runs away exponentially fast in time. 

It is interesting to compare this non-relativistic results with the full relativistic geometry in the presence of the 
cosmological constant. The exact solution of the Einstein equations (the full solution, not just the nonrelativistic 
limit) for the energy-momentum ea. (|2.1l) with p = is the de Sitter universe, 

ds 2 = dt 2 -e 2Ht (dx) 2 . 

Comparing the non-relativistic trajectories ea. (|2.9|l with the exact relativistic geodesies in de Sitter space-time (see 
for example Q) one sees that both exhibit the same exponential runaway behaviour. Therefore, the nonrelativistic 
approximation keeps the essential features of the particle motion in de Sitter space-time. 



III. THE SELF-GRAVITATING GAS WITH A > 



We first consider the hydrostatic description of the self-gravitating gas in thermal equilibrium in the presence of 
the cosmological constant. Next, we present the statistical mechanics of the self-gravitating gas in the presence of the 
cosmological constant and derive the mean field approach to it. The two approaches turn to be equivalent when the 
number of particles tends to infinity. 
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A. Hydrostatics of the self-gravitating gas with A/0 

Let us consider the hydrostatic equilibrium condition [t| for a self-gravitating gas in the presence of the cosmological 
constant. 

VP(q) = -p(q) VV(q) , (3.1) 
where P(q) stands for the pressure at the point q. Let us assume the ideal gas equation of state in local form (see 

MS) 

P(q) = ~ P(q) ■ (3.2) 
m 

Combining eas. If) and (|3.2|) yields the Boltzmann law for the particle density 

p®=poe-9 v ® (3.3) 
where p is a constant, V being the volume of the system. The normalization of density requires that 

d 3 q p(q) = mN . (3.4) 

Inserting eq. l|3.3|l into the Poisson equation 12. 3|) yields, 

VV = 4tt G (po e~9 v ® - 2 A) . (3.5) 



It is convenient to introduce dimensionless coordinates r = q/L where L = V 1 ^ 3 stands for the linear size of the 
system. Ea. ()3.5[) takes thus the form, 



V 2 $(r) + 4tt (rj e* (?) - £) = , (3.6) 



where $(f) = - f V(q) + log %v, 



Cm 2 N T 2 
r)=^ fr - and £ = 2AGm— . (3.7) 

The mass density now becomes p{f) = e*( r ) and the constraint l|3.4|) takes the form, 

d 3 r e* (?) = 1 . (3.8) 



/ 



'unit volume 

The limit for a large number of particles N is well defined provided we choose^ 

N — > oo , L — > oo , A^O, 

with 

— = fixed , A L 2 = fixed . 
L 

It should be mention that an equation similar to eq. 1)3. 6|) appeared in two space dimensions describing; (multi)-vortices 
in the Ginzburg-Landau or Higgs model in the limit between superconductivity of type I and II [Tj. However, for 
the vortex case one has £ = r\ < 0. The negative sign is associated to the fact that like charges repel each other in 
electrodynamics while masses attract each other in gravity. 

In the limiting case £ = ea. 13.6fl becomes the well known Lane-Emden equation in the absence of cosmological 
constant both in the hydrostatic and statistical mechanics approaches 0,0. 

Notice that A plays the same role as the subtraction of a constant density in the Jeans' swindle [l0| . 
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B. The Partition function of the self-gravitating gas with A ^ 

We present now the statistical mechanics of the self-gravitating gas in the presence of the cosmological constant in 
the canonical ensemble. The Hamiltonian is given by eas. (|2.7ll - <|2.8|) . 
The classical partition function of the gas is then, 

1 f f TT d V; d 3 qi -h 



\ (2^) 3 

It is convenient to introduce the dimensionless coordinates rj = % . The momenta integrals are computed straightfor- 
wardly. Hence, the partition function factorizes as the partition function of a perfect gas times the coordinate integral 

Zcoor • 

V N fmT\ 3N/2 

Z =m{-^) z — (3 - 9) 

where 

N 



J\d z n e'" p+ T«»» , (3.10) 



/=i 

and 

JV 



N ^ \n-rA ^ 1 

l<i<j<N 1=1 



where rj and £ are defined by ea. l|3.7|l . 

The free energy takes then the form, 

F = -TlogZ(T,N,V) = -NT log 
We get for the external pressure of the gas, 



eV /mT\ 3/2 
~N~ \ 2^J 



TlogZ coor , (3.11) 



i'=-(w)r^ +T (w) T ^ z -- (3 - 12) 



Evaluating the volume derivative of Z coor with the help of eas. H3.7J) and (|3.1Q|I yields, 



PL = 1+ <U P >-1<^> (3 . 13) 
NT 3NT v ; 



where 



Em, mi , 47rGA \ 
T+-h and U K = — 

Ea. (|3.13|l is the virial theorem for the gas of self-gravitating particles in the presence of the cosmological constant. 

C. Mean field theory 

In the large N limit the 3A^-uple integral in eq. (|3.1U|I can be approximated by a functional integral over the density 
(see refs.Q)- 

Zcoor = I Dp(.)Abe- NS ^ (3.14) 
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with, 



*W->l--I/£$>«"»«-T 



£ y d 3 x p(x) x 2 



+ y d 3 x p(x) In p(x) +ib[l- J d 3 x p(x)] , (3.15) 

and 6 is the Lagrange multiplier enforcing the normalization of the density: 

y d 3 xp(x) = 1 . (3.16) 

The function Z coor thus becomes a functional integral over all density configurations p(x). Using eas, 13.9|l . (|3.14(l and 
(|3.15l) the partition function of the system becomes 



Z = J Dp{.) db e 



where 

F[p(.)} =F + NT S[p(.)} 
stands for the free energy of the gas with density p, while 



(3.17) 



F 



-NT In 



eV ( mT 



is the free energy of the perfect gas. One recognizes the potential energy (the gravitational energy due to the self- 
interaction of the particles and the gravitational energy due to the interaction of the particles with the dark energy) in 
the first line of eq. (|3.15|l , while the second line contains the entropy and the term fixing the total number of particles 
to be N. 



Since the free energy becomes large in the thermodynamic limit (N 3> 1), the functional integral Z coor in ea. (|3.14[) 
is dominated by the extrema of <5>[p(.)]. Extremizing this quantity with respect to the density p(x) yields the saddle 
point equation 

lnp(x) =a + V [ p(y) + ^x 2 . (3.18) 

J |y-x| 3 

This equation defines the mean field approach. We set ib = a + 1 and 

p(x) = e*( x > . (3.19) 
We recognize in ea. (|3.18|l the gravitational potential 

V(x) = f$(x) - a] . (3.20) 

m 

Inserting ea. H3.20f> into ea. l|3.19|l yields the Boltzmann law for the density 

p(x) = e a e"? 7 ' 1 ' , (3.21) 

containing the energy of a particle in the mean field gravitational potential V. e a plays the role of a normalization 
constant. 

Applying the Laplace operator to the saddle point equation l|3.18|l we find the differential equation 

V 2 $(x) + 4tt (77 e*( x ) - A = . (3.22) 

Ea. (|3.22|l is identical to the hydrostatic equilibrium ea. H3.6fl . Therefore, hydrostatics and mean field are equivalent 
in the N — > 00 limit. 

It must be stressed that in the hydrostatic approach one has to assume some local equation of state. Only a 
perfect gas equation yields ea. H3.22|) . In the mean field approach one can compute the local equation of state and 
demonstrate that this is the one of a perfect gas0,0,Q. 
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D. Spherically symmetric case 

For spherically symmetric configurations the mean field equation l|3.6[) becomes an ordinary non-linear differential 
equation. The various thermodynamic quantities are expressed in terms of their solutions. 




0.2 0.4 0.6 0.8 1 1.2 1.4 eta 1.6 



FIG. 1: The pressure at the boundary of the self-gravitating gas without dark energy f(r/, £ = 0) versus rj. PG is the perfect 
gas point, O is the point where the gas collapses (rj° = 1.510 . . .), C is the branch point where c v diverges (rj C = 1.561 . . .), 
and MC is the collapse point in the microcanonical ensemble (r/ MC = 1.259 . . .). 



We consider here the spherically symmetric case where the mean field equation H3.22|) takes the form 

d 2 $ 2 d$ 



dR 2 R dR 



Air (77 e*^ - A = . (3.23) 



We work in a unit sphere, therefore the radial variable runs in the interval < R < R m ax, Rmax = The 

density of particles p(R) has to be normalized according to ea. (|3.16|l . Integrating ea. (|3.23l) from R = to R = R m ax 
yields, 

X] - i = -Rl iax ^(R max ) (3.24) 

Setting, 

£ R / R 

$(R) = u(x)+ln^ , x^^/i^- , (3.25) 
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FIG. 2: The density at the boundary f(r/, £) versus r\ for f = 0, 1, 2.04 . . . , 3. We insert on each graph the point of homogeneous 
density u — and the point of unstability where the isothermal compressibility diverges. The gas is stable from r\ = till the 
point of unstability 



£ R = — and r] R = — , the saddle-point equation (|3.23[l simplifies as, 



2 di 



a U ( X ) 



-1 = 0. 



(3.26) 



dx 2 x dx 

Let us find the boundary conditions for this equation. In order to have a regular solution at origin we have to impose 



u'(0) = 

We find from eas. (|3.24|) - (|3.25|l for fixed values of rj and £, 



r) R -l; R 



(3.27) 



(3.28) 



Eas. H3.27|) and l|3.28|l provide the boundary conditions for the nonlinear ordinary differential equation l|3.26|l . In 
particular, they impose the dependence of uq = u(0) on r] R and £ R . 

Using eas. <|3.19[) and (|3.25|) we can express the density of particles in terms of the solution of ea. l|3.2(j[l 



p(R) = ^ exp[ u (V3^ U 

■q" 



Rr, 



(3.29) 
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FIG. 3: The logarithm of the contrast C(rj, £) versus r\ for £ = 0, 1, 2.04 . . . , 3. Notice that the contrast is unity for r\ — £. 



Recalling that the local pressure is P(R) = P(R) an d usin g ea. (|3.29(l . the contrast C(r),£) between the pressure 
at the center and at the boundary can be written as, 

Using eas. (|3.19|) and i|3.25[l . the pressure at the boundary is given by 



P[Rn 



NT NT N T f R , 

f(v,0 = -ft- PiRmax) = -ft- ^ exp[ W (V3^)] 



V 



V 



V T] R 



while the pressure at the origin takes the value 



P(0) 



NT Z 1 



We find from eq. 



V if 

for small x the expansion of u(x) in powers of x with the result, 

2 4 

u{x) = u + (1 - e u °) *L + e -o ie uo - 1) ^ + 0(x 6 ) . 



(3.30) 



(3.31) 
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where we imposed ea. (|3.27[) . The value of uq follows by imposing ea. l|3.28|) . This power expansion is well suited near 
the center of the isothermal sphere [see sec IV. D]. 

Solving ea. (|3. 26(1 with the boundary conditions u(Q) = Uo, u'(0) = and ea. (|3.28|) yields the function u(R) for the 
whole range of values of r] R and £ . We did that numerically using a fourth order Runge-Kutta method. We plot in 
fig. ^the density at the boundary versus rj for £ = and in fig. [5] the density at the boundary versus 77 for different 
values of £. We plot in fig. 0the contrast C(rj, £) versus 77 for £ = 0, 1, 2.04 . . . , 3 . As we can see from fig. the 
effect of the dark energy is to decrease the value of the contrast in agreement with the repulsive character of the 
dark energy. In particular, the contrast is unity for 77 = £ . In this case the particle distribution is homogeneous [see 
eq. (|4.1() and discussion there] . 

We compute in what follows the thermodynamic quantities as functions of the physical parameters rj and £. 



E. Free energy 




1 1 1 

The free energy vs. eta for xi=0 
The free energy vs. eta for xi=1 
The free energy vs. eta for xi=3 
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FIG. 4: The free energy per particle divided by T , 



versus rj for selected values of £. 



We start by computing the free energy. Using eas. (|3.15() - (|3.17() we find 



F-Fp 
NT 



(3.32) 
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The potential energy vs. eta for xi=0 
The potential energy vs. eta for xi=1 
The potential energy vs. eta for xi=3 
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FIG. 5: The potential energy per particle divided by T , Ep 1 fy^ versus r/ for selected values of £. 



In order to express the Lagrange multiplier a as functions of the physical parameters r\ and £ we use eq. (I3.18|) where 
we can now integrate over the angles with the result, 



- / dR' R 2 e*( fl ') + / dR' R' e*( R ') 

R Jo Jr 



Rmax 



<E>(i?) = a + 4iri] 

We call f(r] R ,£ R ) the density of particles at the boundary R = R max [see ea. p.3Ufl ]. 



2tt ^ , 



tR , 

5_ e <V^) 



Setting R = R m ax in ea. l|3.33|) yields for the Lagrange multiplier 

a = \nf( V R ,S R )-T, R 



where we used the normalization of the densities eq. (|3.1(il) 
Inserting eq. l|3.35[) into the free energy l|3.32[l . we find 



F-F = 1 
NT 2 



hxf(r,,i)-v R -hr 



2 



Rn 



2 7T / R 2 dR 
lo 



$(i?)-|^ 2 



MR) 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



13 




These integrals are computed in appendix A with the result 

F ^ = 3[i-f( V R ,e)]~v R 



NT 



CR / 4 cR\ tli fRrna* 

+ v( i -5^j- 87r ^y dRR2 *w- 

We plot in fig. 0]the free energy as a function of r\ for fixed £. It exhibits two Riemann sheets. We call rjc (£) the 
position of the branch point. The curves have the same qualitative behaviour as for £ = except that the branch 
point T)c(C) increases with £. The dark energy moves the plot of the free energy to lower values. It is a consequence 
of the negative contribution of the dark energy to the energy [ea. (|2.7|> ] and the free energy [ea. i|3.32[l ]. 

When matter dominates over dark energy we have JVm>7A and therefore according to ea. l|3.7|l rj 3> £• In such 
limit the free energy becomes, 

In /(„« f «) - rf + 3 [1 - f(rf, f *)] ■ 
We recognize here the free energy of a self-gravitating gas with N particles of mass m as it must be (see ref.0). 
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F. Energy and entropy 

We compute here the gravitational energy of the particles. The gravitational potential is the sum of the self-gravity 
of the particles plus their interaction with the dark energy [sec cq. <|2.7[l ]. Equivalently, we find for potential energy 
from eqs.fl23, G3, G3 and fiTZty . 

Ep = \ J P(q) V self (q) d 3 q + J p(q) V dark (q) d 3 q = NT [| - \ J d 3 x $(f ) e*& J d 3 f x 2 e*& 



where we use the dimensionless coordinates x = ^ and a is given by eq. 

We compute these integrals for the spherically symmetric case in appendix A with the result, 

%l = 3 [/(r ? «a-i]-^^fi-?CV^5^/(^.a+ 



NT wu ,s ' J 8 s V 5 V R J 4 77^ 

CR rR m ax 

+11 7T ^ / dR R 2 $(i?) . (3.38) 
V Jo 

We plot the gravitational energy as a function of ?y for fixed values of £ in fig. [S] We see a two-sheeted structure with 
the branch point at r/c(£,)- The dark energy pushes the potential energy to lower values. As discussed before for the 
free energy, it is a consequence of the repulsive character of the dark energy. 

We now compute the entropy as S = [E — F]/T using eas. H3.37|) and l|3.38|l for the free energy and the energy, 
respectively, 

^ = 6[/(^a-l]+^-f^ + ^^-(l + ^^)ln/(^a + 

+ 19 tt^ / dRR 2 <P(R). (3.39) 
V Jo 

We plot the entropy as a function of r\ for fixed values of £ in fig. |B| We see again a two-sheeted structure with the 
branch point at tjc{£,)- 

G. Stability of the gaseous phase 

In the precedent paragraphs we present the mean field theory of the selfgravitating gas in the presence of dark 
energy. Mean field theory gives exactly the physical quantities of the gaseous phase in the thermodynamic limit 

N — > 00 , L — > 00 , A^O, 

with 

N 2 

— = fixed , A L = fixed . 

L 

But it does not tell anything directly about the stability of the gaseous phase. 

Necessary conditions for stability in the canonical ensemble are positive specific heat c v and compressibility kt0, 
We compute below c v per particle and Kj- starting from their definitions, 

T 



N VdT L, ' V \ dP 



The point of instability of the self-gravitating gas in absence of dark energy is 77° = 1.510 . . . (see fig. where the 
compressibility diverges^- At this point the isothermal compressibilities has a simple pole and changes sign. It is 
also the point where the speed of sound at the center of the sphere becomes imaginary. Small density fluctuations 
grow exponentially in time instead of exhibiting oscillatory propagation. Such a behavior clearly leads to instability 
and collapse. 
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FIG. 7: The specific heat per particle at constant volume, c v vs. r\ for fixed £ = 0.1, 1. 



In order to compute the specific heat for nonzero £ we introduce the function 

du 



g(x,u ) = y-Q^J ( x ,u ) ■ 
Using eas. (|3.26|) and l|3.27[l we find that this function g obeys to the second order differential equation 

with the boundary conditions 



ff(0,u ) = 1 



O.r 



(o, uo) = o 



(3.41) 
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FIG. 8: The inverse of the isothermal compressibility vs. r\ for fixed £ = 0, 1 and 2. 



c v is computed in the Appendix B with the result 
cv = 3 ( V R 



^)/-?S in /-l^+SS^ 

8 i] H 2 40 rj H 



R 



2 n R 



— tt^t / di? i? 2 uo) 



n 



duo 
drj R 

57 £ fl 1 



9uo 



, 4 , 19^ 



4 »7* (3£*)f 7o 



dx £ 2 <jf(:E, Mq) 



(3.42) 



We see in figHthat the specific heat at constant volume versus r\ for a fixed £ is positive in the first branch till the 
branch point (77 = r\ (£)), where the specific heat tends to infinite and changes its sign. c v is negative in the second 
branch, therefore the second branch of the gaseous phase can not be stable in the canonical ensemble. 



Using eas. <|3.3L)[) the dimensionless isothermal compressibility takes the form, 



Kt 



NT 
~V~ 



K T = 



J {'I > ? I I 3 d^TT 3 gpt 



(3.43) 



3 3 dt 7 

We plot in figlHl versus r\ for fixed £. The isothermal compressibility tends to infinite and becomes negative at the 
point 77 = 77° indicating the collapse of the gas. This unstability happens for values of ry°(£) smaller than 77 (£) (see 
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Table 1). Thus the positivity of the isothermal compressibility is a stronger condition of stability than the positivity 
of the specific heat at constant volume. For £ = Monte Carlo simulations shows that the gas collapse in a very 
dense phase at the point of divergence of the isothermal compressibility ■ Mean field indicates that the gas in the 
presence of dark energy is stable from 77 = till the point 77 = rf (£) . 



a 









1.510... 


1.561... 


1 


1.76... 


1.938... 


1.7 


1.88... 


2.282... 


2.04... 


2.04... 


2.476... 


3 


2.58... 


3.133... 


4.171... 


3.46... 


4.171... 


6 


4.92... 


6.307... 



TABLE 1. 77° (£) and r] C {£_) for different values off 

IV. PHYSICAL PICTURE 

Let us first recall the properties of the self-gravitating gas in absence of dark energy (f = 0) 0] . The self-gravitating 
system can be in one of two phases: gaseous or highly condensed. The pressure at the boundary f(j], f = 0) as a 
function of 77 is given in fig 2] in units of The function f(rj,£ = 0) is calculated from mean field theory and 

describes the gaseous phase, /(jy, f = 0) has two Riemann sheets in the 77-plane. The mean field theory gives exactly 
the physical quantities in the thermodynamic limit. The stability of the gas phase follows by studying the speed 
of sound inside the sphere and was moreover investigated by Monte Carlo simulations 0, 0- The gaseous phase 
is stable only in one piece of the curve f(r), f = 0). This zone of stability depends of the statistical ensemble. In 
the first branch (PG-C) /(t7,0) monotonically decreases for increasing 77 as shown in fig. ^ That is, for decreasing 
temperature or increasing density N/L. The specific heat at constant volume is positive and the gaseous phase can 
be stable in the microcanonical and in the canonical ensemble. In fact, the gaseous phase is stable in the canonical 
ensemble from the point PG (77 = 0, perfect gas) till the point O (77° = 1.510...) where the speed of sound at the 
center becomes imaginary and the isothermal compressibility tends to infinite and changes its sign |2j, |3j . Monte-Carlo 
simulations show that the gas collapses at this point to a extremely condensed phase Q- This is before the point C 
(?7 C = 1.561 . . .) where /(t?, 0) exhibits a vertical slope. C is the branch point of the function /(77JO). The determinant 
of small fluctuations vanishes and the specific heat at constant volume becomes infinite at C|3, 01 • In the second 
branch (C-MC) the specific heat at constant volume is negative. Then, the gaseous phase cannot be stable in the 
canonical ensemble in this second branch but it can be stable in the microcanonical ensemble in (C-MC). In the 
microcanonical ensemble the gaseous phase is stable from PG to MC and collapses in a condensed phase at the point 
MC (t7 MC = 1.259 ...) @,I3. Beyond the point MC /(17, 0) develops a nice spiral form, but the gas is unstable in all 
three statistical ensembles. 

That is, solving just the mean field equation does not tell us whether the gas is stable. One has to compute in 
addition the isothermal compressibility (or the local speed of sound) and the specific heat at constant volume. 

In the gaseous phase the free energy has a minimum for the density p(x), solution of the saddle point equation 
(|3.18l) . This density is the most probable distribution which becomes absolutely certain in the infinite N limit. All 
thermodynamic quantities follow from this density p{x). In the condensed phase, the system becomes extremely dense 
with all particles on the top of each other. The solution of the saddle point equation ceases to be a minimum of the 
free energy and the mean field approach then fails to describe the condensed phase. The condensed phase has been 
found by Monte-Carlo methods in ref.[2|. 

Let us discuss the properties of the self-gravitating gas in the presence of dark energy. 

A. The Triple Point 

The effects of self-gravitation and dark energy are opposed. Self-gravitating forces are attractive while dark energy 
produces repulsion. If the parameters 77 and f , which determine the strength of the gravitation and dark energy, 
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rho(R) vs. R for xi=1 and eta=0.67... (u0=-0.6) 
rho rho(R) vs. R for xi=1 and eta=1 (u0=0) 

rho(R) vs. R for xi=1 and eta=1 .41 ... (u0=0.7) 

1.4 - 



1.3 - 



1.2 - 



1.1 - 



1 - 



0.9 - 




0.8 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 R 0.5 0.6 Rmax .7 

FIG. 9: The densities p(R) vs. the radial coordinate R for f = 1 and r\ = 0.67 . . . , 1 and 1.41 . . .. 

respectively, are equal, a exactly homogeneous sphere p(R) = 1 is a solution (see below). In such special case the 
self- gravitating gas behaves as a perfect gas (with PV = NT everywhere). 

Such homogeneous solution is present for all r\ = £ > 0. There is a peculiar point where the homogeneous solution 
coincides with the point of unstability r\ = rf(£). The point r\ — rf(£) is defined as the point where the isothermal 
compressibility diverges and changes its sign. That is, the two conditions 

/(ry = £,0 = l and V °(0=t 
define an unique 'triple point' 77 = £ = £° = 2.04 . . .. 

B. The critical point as a function of £ 

In fig. [21 we plot the pressure at the boundary versus ry for different values of £. We see that the plots for non zero £ 
exhibit a form analogous to the plot for £ = 0. The effect of the dark energy is to conformally translate the curve to 
higher pressure and higher r\. Notice that the value of the point of unstability ?7 (£) increases for increasing £. For 
£ > the gas collapses at a higher density ^ for a given temperature and for a lower temperature for a given density 
^ than for £ = 0. That is, the domain of stability of the gas increases for increasing £. The dark energy has an 
antigravity effect that disfavours the collapse pushing the particles towards the boundary of the sphere. 
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rho(R) vs. R for xi=2.04... and eta=1 (u0=-1 .41 . 
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FIG. 10: The densities p(R) vs. the radial coordinate R for £ = 2.04. . .. The density increases with R in the first plot. The 
density is uniform in the second plot and coincides with the point of instability r\ = 7?°(£) . 



C. The particle density p(R) 

We study the density p(R) as a function of the radial coordinate R. The density p(R) of the self-gravitating gas 
in absence of dark energy (£ = 0) is always a decreasing function of R. This follows from the attractive character of 
gravitation. At the point PG(?7 = 0) [see fig. the gas is ideal and homogeneous, and p(R) = 1. The farther we 
move from the point PG(?7 = 0) in the £ = case, the more inhomogeneous and faster decreasing is the density p(R). 

The i?-dependence of the density p(R) of a self-gravitating gas in presence of dark energy is more involved be- 
cause the dark energy opposes to the attraction by gravity. The functions p(R) and u(x) have the same qualitative 
dependence in R since they are related by exponentiation [see eq. I|3.29fl ]. 

The behaviour at the center of the sphere (R = 0) is governed by the sign of uo since tt'(O) = [eq. (|3.27[l ] and we 
find from ea. (|33T|) that it"(Q) = |(1 - e u °). Hence, 



signK(O)] 



-sign[u J 



The behaviour at the boundary of the sphere (R = R m ax) is governed by the sign of 77 — £ since according to ea. (13.28(1 

sign[M'(i? maa; )] = -sign[ri R - £_ R ] . 



The particle density exhibits for the stable solutions (namely, for positive kt) one of the following behaviours: 
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FIG. 11: The density p(R) vs. the radial coordinate R is an increasing function for f = 3. 

• decreasing: uq > and r\ > £. The density p(R) decreases from the center of the sphere till the boundary. 
This is the only behaviour present when £ = 0. 

• increasing: uq < and 77 < £. The density increases from the center of the sphere till the boundary. 
More precisely, we find for the stable solutions the following properties in the interval < x < y3^: 

• when 77 > £, is positive and u'(x) < 0. 

• when 77 < £, u(ar) is negative and u'(a;) > 0. 



The case 77 = £ is particularly interesting, because we have u'(0) = and u' ^v^ fl J = 0- This implies that 

u(x)=0, (4.1) 

is a solution of ea. (|3.2(il) . However, there exist another solution u ^ for r\ = £ > 1.638 . . ., but this solution is never 
stable since kt < for it. For the u = solution, p(R) = 1 and the gas is homogeneous. The self-gravity of particles 
is exactly compensated by their interaction with the dark energy. We notice that the u = homogeneous solution is 
stable only for < £ < 2.04 . . .. 

We recall that in the vortex case 77 = £ < there are nontrivial solutions for 7j(i?)^lJ- 
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Let us now study the particle density p(R) as a function of £. There are three stable cases (kt > 0). We illustrate 
these cases plotting f(r), £), the surface density in units of in fig. That is, f(r), £) vs. 77 for £ = 1, 2.04 . . ., and 
3. The densities p(R) are plotted as functions of R for £ = 1, 2.04 . . ., and 3 and chosen values of 77 in figs. iQltolTTl 

• First case < £ < 2.04 . . . 

We illustrate this case choosing the value £ = 1. As we see in fig. [21 the line 77 = 1 intercepts ,/(?/,£ = 1) 
just in one point. This intersection corresponds to the solution u(x) = 0. At this point £ = 77 we also have 
mo = 0. We want to notice that when we move over the curve f(r], £ = fixed) increasing 77, uq always increase. 
Therefore, uq < above the point where u(x) = while uo > below the point where u(x) — 0. This property 
holds for all values of £ (for £ > 1.638 . . . the line 77 = £ intercepts f(j], £) vs. 77 in more than two points, for 
1.638 ...<£< 2.04 . . . the point u = is in the upper branch). 

We have for this case r/ (£) > £ . There are two zones of stability in fig. 2 for f(r), £ = 1): 

-the zone between the point 77 = and the point of homogeneous density (77 = £, u — 0), where uq < and 
77 < £ (< 77 (£)), which corresponds to a density increasing with R. 

-the zone between the point of homogeneous density (77 = £, u — 0) and the point of unstability 77 = 77°(£), where 
uo > and £ < 77 < ?7 (£), which corresponds to a density decreasing with R. 

The gas is unstable in the zone beyond the critical point of unstability 77 = 77 (£) . 

We plot in fig. OH the densities p(R) vs. R for £ = 1 and 77 = 0.67 . . . , 77 = 1 and 77 = 1.41 . . .. The first plot 
corresponds to a increasing density. The dark energy dominates over gravitational attraction. The second plot 
corresponds to the uniform density associated to u{x) = 0. The density is homogeneous and dark energy exactly 
compensates gravitational attraction. The third plot corresponds to a decreasing density. The gravitational 
attraction dominates over dark energy. 

• Second case £ = £° = 2.04 . . . 

This is the value of £ where rf{£) — £. This is a 'triple point' since it is defined by two conditions 

/fa = £,0 = 1 and 77°(£)=£ 

In this case the point of unstability 77 = ?7 (£) coincides with the point of homogeneous density (7/ = £°, u = 0). 
In fig. 2 /(?7, £ = 2.04 . . .) the gas is stable in the zone between the point 77 = and the point of homogeneous 
density (77 = 77°(£°) = £ , u = 0), where uq < and 77 < 77 (£ ) = £ , which corresponds to an increasing density. 
The gas is unstable in the zone beyond the critical point of unstability 77 = 77 (£). 

We plot in fig. ^|the densities p(R) vs. R for £ = 2.04... and 77 = 1, r] = 2.04 . . .. The first plot exhibits a density 
increasing with R. The dark energy dominates over the gravitational attraction. The second plot corresponds 
to the point of homogeneous density u(x) = 0, which coincides with the point of unstability 77 = 77 (£). 

• Third case £ > 2.04 . . . 

We illustrate this case choosing the value £ = 3. We have ?7 (£) < £ for this case. The homogeneous solution 
u = is unstable because kt < for it. In the graph /(?7, £ = 3) the gas is stable in the zone between the 
point rj = and the point of unstability 77 = 77 (£), where uq < and 77 < 77°(£)(< £), which corresponds to an 
increasing density. 

We plot in fig. ^]the densities p(R) vs. R for £ = 3 and 77 = 1.379 . . .. The plot corresponds to an increasing 
density. The dark energy dominates over gravitational attraction. 

It could be mentioned that the homogeneous solution u(x) = that we have for 77 = £ is the nonrelativistic analogue 
of the static Einstein universe. It was precisely thanks to the introduction of the cosmological constant that Einstein 
found such static and homogeneous universe. 



D. The Diluted expansion 



We can obtain power expansions in 77 and £ for the various physical quantities from eqs. 
analytic expressions are valid for 77 <§C 1 and £ <^ 1. We find, 



V 



3 / ^ 51 ^2 I 83 
H (77 — £) H £ 2 - 

10 w v 1400 ^ 



33 

77 £ + rj 2 + Of?? 3 , 77 2 £, 77 £ 2 , £ 3 ) 

1400 350 ' W ' ' 



and ipOT)) . These 
(4.2) 
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Using now eas. (|3.30|) . (|3.31(l and i|4.2[l . we obtain 

f(v,o = i - ! ^-i^( 7 ?+ 2 o +o(vWz,v e,e) (4.3) 

Ea. (|4.3() reduces for £ = to the corresponding expression in absence of dark energy (see |2J). For rj = £ we retrieve 
the solution u — with / = 1. We see comparing with figEl that the first order approximation f(r),£) = 1 — 2^ is 
valid even for large r\ and £ till the vicinity of the point where the slope becomes vertical. This unexpectedly large 
domain of validity may be explained by the smallness (y^g) of the second order coefficient. 

We clearly see that the cosmological constant translates as a whole the plot of f(r], £ = 0) by an amount +£. That 

is, 

for rj < r]c(0- 



V. THE LIMITING CASE £ » 77 



We consider here the limit where the dark matter dominates over the ordinary matter (£ ^> 77). We have a gas of 
particles where their self-gravitational interactions are dominated by the interaction with the dark energy. 
The mean field equation (|3.22|) becomes linear 

V 2 $(x) = 4tt£ . 



The exact solution for spherical symmetry takes the form, 



R 



Rmax 



(5.1) 



where a is independent of R. We recall that R runs in < R < R ma x with R n 
The density of matter follows from eas. I|3. 19(1 and (|5.1(1 . 

p(R) = e a eVlssi) 

Notice that this density increases monotonically with R [see sec. IV. C]. 
The contrast is given by 



P(0) 



p(Rn 



= p 2 



The density normalization l|3.16|l imposes that 



dy y r e 2 y 



This integral can be expressed in terms of the probability integral $(2;) as 

follows|I3 



-<*(€) 



where 



<f>(x) 



dy e" 



(5.2) 



(5.3) 



It is instructive to derive these results dropping the particle self-interaction in the partition function (|3.9() . The 
Hamiltonian is in this limit. 



H = J2 



2mj 
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The coordinate partition function Z coor then factorizes and we find from eq. 

rN / m rr\ 3N/2 



V N fmT\ 



where a(£, R ) is given by ea. (|5.3|) after using the dimensionless coordinates r = q/L. 
The density at the boundary (|3. 34|) takes here the form 

/((),£*) - p{R max ) = 

We plot in fig. 1121 the density at the boundary versus £. We give below the thermodynamic quantities in terms of 
/(0,£ ) and £ R . That is, free energy, potential energy, pressure at the boundary, entropy, specific heat at constant 
volume and isothermal compressibility of the gas: 




xi A R 16 



FIG. 12: The density at the boundary f(rj R = 0,£ R ) versus £ R when the effects of the dark matter dominate over the 
selfgravitation. 
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= \ + 1 e m t*) + § [i - m ^ [i + §/(o, n 



i 



/(<U fl ) 



(5.4) 



i^o and So are respectively the free energy and the entropy of the perfect gas. The second and third relations in 
eas. (|5.4t yield virial theorem, 



PV 



1-- 3 E P 



which is the limiting case of ea. l|3.13|l for 77 = 0. 

We checked that c v and kt are positive for all £ R > 0. The gaseous phase is thus stable in this limit. 

We compute the asymptotic behavior of the thermodynamic functions for £ R — ► 00 
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(5.5) 



and for 



E L = _3{f_3^ + 0(£S . ) 
NT 10 175 ^ y 



o(^ J ) 
s-s 

N 



- 2 ' 175 + «*> ' 
The particle density for £ R — > 00 takes the form, 



7VT ~~ 10 350 



3^ 



350 



KT = 1 - g + + 0(^ 3 ) 

15 1575 



p(i?) 



_^±1 -^fi-(^) 2 ] 



(5.6) 



Notice that the £ fl — > 00 corresponds to keeping A (and if) fixed in the thermodynamic limit [see ea. l|l.l|) ]. 



VI. DISCUSSION AND CONCLUSIONS 



The behaviour of the self-gravitating gas is significantly influenced by the cosmological constant (or not) depending 
on the value of the ratio 

R A = 1 = 2AV = 2 ( dark energy \ = ^ A_ 
r) m N \ mass J v py 

where py = N m/V is the average mass density in the volume V and where the subscript Q v indicates mass and 
dark energy inside the volume V. 

We display in Table 2 the typical value of Ra for some relevant astrophysical objects. Except for the universe as 
a whole Ra takes small values. For such values of R\ the effect of the cosmological constant is negligible as we see 
from the previous sections. The smallncss of R\ stems from the tiny value of A = 0.663 10 _29 g/cm 3 [13j. The particle 
mass density is much larger than A in all situations except for the universe as a whole. However, the nonrelativistic 
and equilibrium treatment does not apply for the universe as a whole. 
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Object 


Ra 


Stars 


~ 1(T 29 


Supergiants 


99 

~ 10~ 22 


Galaxies 


- 1(T 6 - 1(T 5 


Clusters of Galaxies 


~ lO" 2 


The Universe 


~ 4 



TABLE 2. The ratio i? A = £ for some astrophysical objects 

For the cold clouds in the interstellar medium (ISM) considered in refs.0,|]JQ it is more convenient to express £ as 

£=1.210- 6 ^[L(rf 

where the mass of the particles m is expressed in units of the hydrogen mass and T in Kelvins. For typical cold clouds 
in the ISM we have 1 < L(pc) < 100 and T ~ 5 ; 3]. This yields 



£~10~ 3 to 10~ 7 



while r\ < 1. Again, we find Ra <C 1. 



The isothermal sphere in the presence of a cosmological constant has been considered in refs.pj). The authors in 
ref.^4[ restrict themselves to the case where the gravitational force points inward in all points. This happens when 
r\ > £ (using our notation). As showed here, when the dark energy dominates, there are stable solutions for £ > 77 
where p(R) increases with R. In these situations where the dark energy dominates the gravitational force points 
outward in all points. 

It must be noticed that the mean field equation (|3.26|) admits solutions for which the density is not a monotonous 
function of R. However, such solutions are unstable since they correspond to negative values for the compressibility 

KT- 

The effect produced by the cosmological constant is of an opposite nature to the Jeans' or gravothermal instabilities 
where the self-gravitating gas collapses due to the mutual attraction. Namely, the points C and MC in figs. 1 and 2. 



VII. APPENDIX 
A. Appendix A 

The goal of this appendix is to compute in the spherically symmetric case the following integrals 

2 7T 



A = J d 3 xe^ 
Using ea. H3.6|) . we obtain 
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The terms in the first lines of A and B can be integrated twice by parts with the following result in the spherically 
symmetric case, 



1 = c-^^^-tO-I^ 



I + 2J 



B 



I - — 1 lnf( V ,0- — fl-^O 

2 77/ JW ' W 4 V 5? ? V 2 



where f(rj,£) is defined by ca. (|3.34() and 



di? i? 2 



(7.1) 



(7.2) 



V JO '/ JO 

We compute now / in terms of J, /, 77 and £. Let us consider the function 

B(R) = R 3 e* (i?) . 
We take the derivative with respect to R and use ea. (|3.23[l . 

B'(R) = 3 i? 2 e* (i?) - -i- i? 3 $'(i?) -J—R 2 \^'{R)} 2 + R 3 - $'(JZ) . 

47r?7 27T77 77 

We now integrate B'(R) between and R m ax- Integrating by parts the second and the fourth terms and using 
eas.lj3.8JI and (|7.2|> yields, 



I = -6 J + 6 [1 - f{r,, 0] + 6 1 In /fa, £) - tj* + (2 - ^ 

77 \ V 



(7.3) 



We finally get for A and £? using eas. (|7.1|) and (|7.3|) 



-1 = «[ 1 -./(</■ a] + (1 + 40^/(^0-^ + ^0-^) - 1/ 



5 = 6[/(t,,0-1]- 1 



ill 

2r? 



4 10 77/ 2 



(7.4) 



B. Appendix B 

In this appendix, we will compute the specific heat at constant volume. Using eas. l|3.7|l . (|3.39|) and (|3.4U|I . we find 



tR g 
+207r ^ T df 



V 



dR R 2 $(R) 



v Jo 

Using eas. (|3.3U|) we calculate the partial derivative of / with respect to T, 



r 



^ -Tl 



<9r 



We obtain from eqs.|(23, (jiHgl) and 



v 



t( 9 A 

\9T) V 



Using eqs 



. CTTt . IT2Tt and (SH|| we 



find 



<9T 



v Jo 



(7.5) 



(7.6) 
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+ JL t ( — ^ -J: 



/3 



da; a: 2 g(x, uq) 



Integrating by parts the first integral of the right member of the last equation, we obtain 



dT 



dR R 2 &(R) 



v Jo 



/3 £ R 



_JL„,^.„ 0) + JL r ( 

Using eas. (|3.25|) . I|7.6|l and l|7.8|l we express c v as in eq 



9 1 

dup \ 1 

dT) v (3C«) 3 / 2 7o 



dx X 2 u(x, Uq) 

dx x 2 g(x, Uq) 



(7.7) 



(7.8) 
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